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ABSTRACT 

Direct numerical simulations of turbulent Hall dynamos are presented. The evolution of an initially 
weak and small scale magnetic field in a system maintained in a stationary turbulent regime by a stirring 
force at a macroscopic scale is studied to explore the conditions for exponential growth of the magnetic 
energy. Scaling of the dynamo efficiency with the Reynolds numbers is studied, and the resulting total 
energy spectra are found to be compatible with a Kolmogorov type law. A faster growth of large scale 
magnetic fields is observed at intermediate intensities of the Hall effect. 

Subject headings: MHD — magnetic fields — stars: magnetic fields — stars: neutron — accretion disks 



1. INTRODUCTION 

In recent years, the relevance of two-fluid effects 
has been pointed out in several astrophysical (Bal- 
bus & Terquem 2001; Sano & Stone 2002; Mininni, 
Gomez, & Mahajan 2002, 2003a) as well as labora- 
tory plasmas (Mirnov, Hegna, & Prager 2003). The 
standard magnetohydrodynamic (MHD) framework 
for the study of astrophysical plasmas may not be ade- 
quate in the presence of strong magnetic fields and/or 
low ionization; the electric conductivity then is not 
isotropic and nonlinear effects arise in Ohm's law. In 
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low temperature accretion disks around young stellar 
objects or in dwarf nova systems in quiescence, for 
example, the plasma is only partially ionized with a 
small abundance of charged particles. As a result, two 
new effects appear in Ohms's law: ambipolar diffusion 
and Hall currents. The predominance of either of these 
effects is determined by the ionization fraction and the 
plasma density (Sano & Stone 2002; Braginskii 1965). 
The impact of ambipolar diffusion on dynamo action 
has been already studied by Zweibel (1988) (see also 
Brandenburg & Subramanian (2000)). The Hall effect 
affects the dynamics of protostellar disks (Balbus & 
Terquem 2001). In neutron stars magnetic fields are so 
strong that the Hall term can be even more important 
than the induction term for the magnetic field evolu- 
tion (Muslimov 1994). The Hall effect is also known 
to be relevant in others astrophysical scenarios (Sano 
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& Stone 2002; Mininni, Gomez, & Mahajan 2003b, 
and references therein). 

The Hall effect, when strong, is expected to seri- 
ously affect the MHD results on the generation of mag- 
netic fields in astrophysical and laboratory plasmas by 
inductive motions in a conducting fluid (dynamo ef- 
fect). More specifically, it is expected to modify the 
growth and evolution of magnetic energy, since the ad- 
dition of the Hall term to the MHD equations leads 
to the freezing of the magnetic field to the electron 
flow (in the non-dissipative Umit) rather than to the 
bulk velocity field. The first studies on the impact of 
Hall currents on dynamo action (Helmis 1968; Helmis 
1971) were carried out using mean field theory and the 
first-order smoothing approximation (Krause & Radler 
1980). Helmis obtained decreasing dynamo action as 
the strength of the Hall terms increased. Recently, 
the impact of the HaU effect has been studied using 
the kinematic approximation or focusing on particu- 
lar geometries, both in numerical simulations (Galanti, 
Kleeorin, & Rogachevskii 1994) and in experimen- 
tal and theoretical studies (Heintzmann 1983; Ji 1999; 
Rheinhardt & Geppert 2002; Mirnov, Hegna, & Prager 
2003). A general closure scheme was also proposed 
to compute the contribution of the Hall term to the 
dynamo action a-effect (Mininni, Gomez, & Mahajan 
2002) using mean field theory and the reduced smooth- 
ing approximation (Blackman & Field 1999); it was 
found that the Hall effect could either suppress or en- 
hance dynamo action. Expressions for the turbulent 
diffusivity using this closure were also derived for par- 
ticular cases in Mininni, Gomez, & Mahajan (2003a). 

In the present work, we report results from di- 
rect numerical simulations of the dynamo action in 
MHD and Hall-MHD at moderate Reynolds num- 
bers with strong kinetic helical forcing. The study 
of the scaling of the dynamo efficiency with increas- 
ing Reynolds numbers is the main aim of this pa- 
per. In a previous paper (Mininni, Gomez, & Ma- 
hajan 2003b), we showed that three distinct dynamo 
regimes can be clearly identified: (1) dynamo ac- 
tivity is enhanced (Hall-enhanced regime), (2) it is 
inhibited (Hall-suppressed regime), (3) it asymptot- 
ically approaches the MHD value (MHD regime). 
These regimes arise as a result of the relative order- 
ing between the relevant lengthscales of the problem. 
Namely, the energy-containing scale of the flow, the 
Hall length, and the correlation length of the mag- 
netic seed. Simulations in Mininni, Gomez, & Ma- 
hajan (2003b) were performed with 64^ spatial grid 



points and fixed Reynolds numbers. This study left 
unanswered the question whether the enhancement of 
dynamo action by HaU effect would increase or de- 
crease with increasing Reynolds numbers and scale 
separation. 

Theoretical estimates using mean field theory and 
a particular choice for the small scale fields suggest 
that the efficiency of Hall-MHD dynamos compared 
with the MHD counterpart will increase as the scale 
separation is increased (Mininni, Gomez, & Mahajan 
2002). In this work we present direct simulations with 
higher spatial resolutions and Reynold numbers that 
confirm this result. Also, the magnetic, kinetic, and 
total energy spectra developed in Hall-MHD turbu- 
lence are calculated and studied. These spectra corre- 
spond to the dynamo regime, where no imposed cur- 
rents are present, i.e, the magnetic fields are purely 
self-generated. This regime is rather general and prob- 
ably common to most astrophysical flows (Haugen, 
Brandenburg & Dobler 2004). 

The paper is organized as follows. In Section 2 
we present the general equations describing the evo- 
lution of the fields. In Section 3, the code used to nu- 
merically integrate the Hall-MHD system is described, 
and all the simulations made with different resolu- 
tions and parameters are listed. Section 4 presents the 
results obtained in the MHD limit of our equations; 
the results obtained are similar to those of other au- 
thors (Meneguzzi, Frisch, & Pouquet 1981; Branden- 
burg 2001), and provide a reference set to compare 
the Hall-MHD solutions with. Section 5 is devoted to 
the results of Hall-MHD simulations. In Section 6 we 
discuss a subset of MHD and Hall-MHD simulations 
made to study the evolution of large scale magnetic 
fields. Finally, in Section 7, we give a brief summary 
of the current effort. 

2. THE HALL-MHD SYSTEM 

Incompressible Hall-MHD is described by the mod- 
ified induction and the dissipative Navier-Stokes equa- 
tion, 

— = Vx [(I7-eVxB) xB] -hr?V2B(l) 
^ = -{U-V)U+{B.V)B- 

-V (^P + ^^ + F + uV^U (2) 
V • B = = V • C/ , (3) 
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where F denotes a solenoidal external force. The ve- 
locity U and the magnetic field B are expressed in 
units of a characteristic speed Uo, e measures the rela- 
tive strength of the Hall effect, and 77 and i' are the (di- 
mensionless) magnetic diffusivity and kinematic vis- 
cosity, respectively. Note that the measure of Hall ef- 
fect e can be written as 



e = 



LhoM 



(4) 



where Lq is a characteristic length scale (the size of 
the box in our simulations is equal to 27r), and the Hall 
length 

LHall = TT" ) (5) 

LOpi Uo 

is given in terms of the characteristic speed Uq and the 
characteristic Alfvenic speed t//i. In particular, we are 
free to choose Uo = Ua as our characteristic velocity, 
reducing Lnaii to the ion skin depth. 

Equation (5) is valid in a fully ionized plasma. We 
will work under this assumption without any loss of 
generality. In the more general case of partially ion- 
ized plasmas, the values of Lnaii and e are those re- 
ported by Sano & Stone (2002), or in Mininni, Gomez, 
& Mahajan (2003b). Typical values of e in astro- 
physics are also mentioned in these papers. 

To keep in mind astrophysical scenarios, we just 
recall three examples. In a protostellar disk the Hall 
scale is larger than the dissipation scale typically by 
two orders of magnitude, but smaller than the largest 
scales of the system (Balbus & Terquem 2001). There- 
fore we expect e < 1 but with Lnaii larger than the 
dissipation scale. In some dwarf nova disks and proto- 
planetary disks e « 1 (Sano & Stone 2002). As pre- 
viously mentioned, in neutron stars e > 1 (Muslimov 
1994). 

The Hall-MHD system has three well-known ideal 
(n = v = 0) quadratic invariants 



E 



Hrr 



A-BdV ^ 



(6) 
(7) 



K 



^ j{B + euj) ■ {A + eU)dV . (8) 



Here E is the energy, Hm is the magnetic helicity, and 
K is the hybrid helicity, which replaces the cross heUc- 
ity from magnetohydrodynamics. The vector potential 
A is defined by B = V x A, and w = V x [/ is the 



vorticity. Conservation of these ideal invariants during 
the evolution of the system provides a check on the 
simulation. 

3. THE CODE 

The pseudospectral code used in Mininni, Gomez, 
& Mahajan (2003b) was modified to run in a Beowulf 
cluster using MPl. We integrated the Hall-MHD equa- 
tions (l)-(3) in a cubic box with periodic boundary 
conditions. The equations were evolved in time using a 
second order Runge-Kutta method. The total pressure 
Pt = P + /2 was computed in a self-consistent 
fashion at each time step to ensure the incompressibil- 
ity condition V C/ = (Canuto, Hussaini, Quarteroni, 
& Zang 1998). In Fourier space, taking the divergence 
of equation (2) we obtain 



PT(fe) = ^fe • 



(C7 • VC/)fe - {B ■ VB)fe 



(9) 



where the hat denotes a spatial Fourier transform, and 
k is the wavenumber vector. 

To satisfy the divergence-free condition for the 
magnetic field, the induction equation (1) was replaced 
by an equation for the vector potential 



eV X B) xB + eVpe 



■r]V^A 



(10) 



where p,, (electron pressure) was computed at each 
time step to satisfy the Coulomb gauge V • A = 0, 
solving an equation similar to equation (9). 

We present results from different runs with 77 = 
0.05,77 = 0.02, and 7? = 0.011. For the first value of 77, 
simulations with 64^ and 128^ spatial grid points were 
performed to check convergence. The rest of the sim- 
ulations were made with 128'^ grid points (77 = 0.02), 
and 256^ grid points (rj — 0.011). All the runs were 
made with magnetic Prandtl number ly/rj = 1. There- 
fore, hereafter we will only consider a single Reynolds 
number (i.e. both kinetic and magnetic), defined as 



R = 



ULq 
V 



In the study of turbulent flows, the number 

V 



(11) 



(12) 



constructed from Taylor's length scale [A = {{U^) / (w^))!/^] 
is often considered. Note that this definition of 
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Taylor's micro-scale might differ from other defini- 
tions (for instance, in connection with experiments 
on fluid turbulence) in factors of order unity (Pope 
2000). These two Reynolds numbers were respec- 
tively w 100 and R\ « 20 for the first value of ly 
and rj, RKi 300 and R\ w 40 for the second one, and 
R « 560 and Rx « 60 in the last case. The energy 
injection rate was approximately the same for all these 
runs. 

The simulation begins by subjecting the Navier- 
Stokes equation to a stationary helical force F (given 
by eigenfunctions of the curl operator) operating at 
a macroscopic scale k force = 3 (Mininni, Gomez, 
& Mahajan 2003b) to reach a hydrodynamic turbu- 
lent steady state. The resulting statistically steady 
state is characterized by a positive kinetic helic- 
ity. The relative heUcity in runs with R = 300 is 
2 iJfc/((C/2) (w2))V2 ~ 0.4, and this value decreases 
slightly for larger Reynolds numbers. The kinetic he- 
licity is defined as 

Hk = \ju-udV. (13) 

Once the hydrodynamic stage of the simulation 
reaches a steady state, a non-heUcal but small magnetic 
seed was introduced. This initial magnetic seed was 
generated by a J-correlated vector potential centered at 
kseed = 13 for the = 100 runs and kgeed = 35 for 
the R = 300 and R = 560 simulations. The run was 
continued with the same external helical force in the 
Navier-Stokes equation, to study the growth of mag- 
netic energy due to dynamo action. 

Another set of simulations was made under the 
same conditions but with v = rj = 0.02 and k force = 
10 (128^ grid points), to study the changes in the 
growth of the large scale magnetic field in the presence 
of the Hall effect. In this case the Reynolds numbers 
are smaller (R w 220 and R\ « 20), and the turbu- 
lence is weaker, since there are not enough modes in 
Fourier space for a direct cascade to develop properly. 
On the other hand, there are more Fourier modes to 
study the inverse cascade and the growth of large scale 
fields. The results of these simulations are discussed 
in Section 6. 

In all our simulations the Kolmogorov's kinetic and 
magnetic dissipation length scales were properly re- 
solved in the computational domain, i.e. we made sure 
that the dissipation wavenumbers remain smaller than 
the maximum wavenumber allowed by the dealiasing 
step, namely kmax = 128/3. 



4. MHD DYNAMOS 

In this section we briefly present the results from 
MHD simulations. The results are in good agree- 
ment with previous simulations of dynamo action un- 
der periodic boundary conditions (Meneguzzi, Frisch, 
& Pouquet 1981; Brandenburg 2001). These simula- 
tions are intended for comparison with the Hall-MHD 
runs, and therefore some specific results of the MHD 
simulations are discussed in more detail in the follow- 
ing sub-sections. 

4.1. Magnetic energy evolution 

In Figure 1, we show the magnetic and kinetic en- 
ergy as a function of time in MHD runs (e = 0) for 
different Reynolds numbers. The turnover time for 
i? = 300 is T = 27r/{kforce (U^y^^) ~ 0.3. 

Two phases can be clearly identified in the evolu- 
tion of the magnetic energy. After a first stage with 
exponential growth (which can be considered as the 
kinematic dynamo stage), the magnetic energy satu- 
rates and reaches equipartition with the kinetic energy 
(see Figure 1). In the first stage, the magnetic energy is 
still weak and the velocity field is not strongly affected 
by the Lorentz force. Note that during this exponential 
growth of magnetic energy, the kinetic energy remains 
approximately constant. 

This kinematic dynamo stage can be understood at 
least qualitatively considering the mean field induction 
equation (Krause & Radler 1980) 

— = Vx (UxB + aB) + rjeffV'^B. (14) 
Here the overline denotes mean field quantities, and 
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Fig. 1. — Kinetic energy Ek (above), and magnetic 
energy Em (below) as a function of time (e = 0). 
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r]eff is the magnetic plus turbulent diffusivity. The 
MHD a-effect (Pouquet, Frisch, & Leorat 1976) 



to confirm this result although in some cases the spec- 
trum tends to be a little shallower. 



a = -{-u-Vxu + b-yxb) , (15) 

represents the back-reaction of the turbulent motions 
in the mean field, and gives exponential growth of 
magnetic energy in the kinematic regime for helical 
turbulence. Here u and b are respectively the fluc- 
tuating velocity and magnetic fields, and r is a typ- 
ical correlation time for the turbulent motions. At- 
tempts to measure this quantity in direct simulations 
were made by Cattaneo & Hughes (1996), and Bran- 
denburg (2001). 

As the Reynolds number increases, the saturation 
field strength increases, although it seems to reach an 
asymptotic value. After the saturation, the magnetic 
energy keeps growing slowly on a resistive timescale 
(Brandenburg 2001). This late growth takes place 
mainly at large scales, as will be shown in the energy 
spectrum. As far as we know, there are no simulations 
of MHD dynamos with periodic boundary conditions 
showing generation of large scale fields on shorter 
times. 

4.2. Energy spectrum 

Figure 2 shows the kinetic and magnetic spectra 
at different times for a run with R = 300. In the 
early stages the magnetic energy grows uniformly at 
all wave numbers. After the saturation (t « 5) the 
emergence of a large-scale field can be clearly seen 
in the spectrum. At t « 18.4, when the system has 
already reached equipartition (see Figure 1), the mag- 
netic energy at large scales (small wave numbers) still 
keeps growing, albeit slowly. As a result, the large 
scale magnetic field reaches super-equipartition with 
the kinetic energy. An excess of magnetic energy can 
be also observed at small scales. 

The slope of the total (magnetic and kinetic) energy 
spectrum in the inertial range is consistent with Kol- 
mogorov's k^^/'^ law and in good agreement with sim- 
ulations of helical MHD turbulence with higher spatial 
resolutions. Kida, Yanase, & Mizushima (1991) found 
that the total energy spectrum in MHD dynamo simu- 
lations is of the form 

E{k) = CKS^'^k-^'^. (16) 

where e is the total dissipation rate, and Ck is a Kol- 
mogorov's constant. Simulations at higher spatial res- 
olutions (Haugen, Brandenburg & Dobler 2003) seem 



4.3. Magnetic helicity 

The magnetic helicity H„i is displayed in Figure 
3. The initial magnetic field is non-hehcal, but during 
the dynamo process net magnetic helicity is generated 
with a sign opposite to that of the kinetic helicity. This 
hehcity is located mostly in large scale structures, as 
will be shown in Section 6. 

From mean field equations we obtain for the large 
scale magnetic helicity (Mininni, Gomez, & Mahajan 
2003b) 

^ = 2 y - 77e// J ■B)dV, (17) 

where J = V x B is the electric current density. 
This equation represents a transfer of magnetic helic- 
ity from small scales to large scales. The mean field 
helicity grows with the same sign as the a coefficient 
(opposite sign as the kinetic helicity). Our results are 
in good agreement with this relation, as well as previ- 
ous simulations of MHD dynamo action (Brandenburg 
2001; Mininni, Gomez, & Mahajan 2003b). 

Although this generation of magnetic helicity by 
dynamo action is expected to decrease as the magnetic 
Reynolds number Rm increases, in simulations with 
higher Reynolds numbers the growth of magnetic he- 
licity seems to reach an asymptotic value. 

5. HALL DYNAMOS 

To quantitatively assess the role of the Hall effect 
on dynamo action, we display results from runs with 
different values of e : the MHD run (e = 0) and Hall- 
MHD runs with e = 0.066, 0.1 and 0.2. We will fo- 
cus on the Hall-enhanced dynamo regime (Mininni, 
Gomez, & Mahajan 2003b). Note that all these val- 
ues of e correspond to dynamos where the Hall ef- 
fect is only relevant in a fraction of the scales in- 
volved. The Hall inverse length scale for these runs 
is measured by knaii = 15, 10 and 5 respectively 
(kHaii = 1/e). All length scales smaller than the Hall 
scale are expected to be strongly affected by the Hall 
effect. The Kolmogorov's kinetic dissipation scale 
[k^ = ((a;2) /iy2)i/4] is 20 when R = 100, 

w 40 when R = 300, and k^ w 75 when R = 560. 
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5.1. Magnetic energy evolution 




Fig. 2. — Mean kinetic energy spectrum (thick line), 
total energy spectrum (thick dashed Une), and mag- 
netic energy spectrum at different times (e = and 
R — 300). The Kolmogorov's slope is shown as a ref- 
erence. 




Figure 4 shows the kinetic and magnetic energy as 
a function of time for the MHD and a Hall-MHD run 
with R = 300 and e = 0.1. At early times, the evo- 
lution of magnetic energy in MHD and Hall-MHD is 
similar. Using mean field theory, the induction equa- 
tion for the mean magnetic field reduces to 



dB 

dt 



Vx [{U -eVxB) xB + aB]+r]effV^B, 

(18) 

with the HaU-MHD a-effect now given by (Mininni, 
Gomez, & Mahajan 2002) 



a = - (-w^ • Vxti« -I- 6 • Vx6 



-e6- VxVxw^) . 



(19) 



Here tx^ = ti — e V X 6 is the fluctuating electron flow 
velocity. When the fluctuating magnetic field is weak, 
this expression reduces to equation (15) and the Hall 
effect can be dropped. Therefore, the first stage corre- 
sponds to a kinematic dynamo during which the HaU 
effect is negligible. 

After this stage, and when the dynamo-generated 
magnetic fields are strong enough for the Hall effect to 
become non-negligible, the evolution changes and the 
magnetic energy keeps growing but at a different pace 
(see also Figure 5). A third stage can be identified, 
when the velocity field is affected by the Lorentz force 
and the dynamo reaches saturation. The kinetic energy 
drop is not as intense as in the MHD case, and the in- 
crease of magnetic energy in this final stage is larger 
than in the MHD case for moderate values of e. Fi- 
naUy, a state with more magnetic energy than its MHD 
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Fig. 3. — Magnetic helicity {—Hm)- 



Fig. 4. — Magnetic (below) and kinetic energy (above) 
as a function of time for two runs with e = 0.1 and 
e = (i? = 300). 
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counterpart is reached (by a factor 2.3 when R = 300 
and e = 0.1). 

Figure 5 shows the evolution for a Hall-MHD run 
with R = 300 and e = 0.2. Here, after the expo- 
nentially growing stage, the magnetic field saturates at 
an ampUtude smaller than the previous run. Note that 
equipartition between the kinetic and magnetic energy 
is not reached. 

When a turbulent stationary state is attained, the 
sum of magnetic and kinetic energy is not equal to the 
initial kinetic energy. This is related to the fact that, 
when the magnetic seed is introduced, a new channel 
for energy dissipation arises, 

^ = -u j uj^dV -ri j J^dV , (20) 

As in previous simulations (Mininni, Gomez, & Maha- 
jan 2003b), it is found that the final energy reached for 
the Hall-MHD runs is larger than the value obtained 
for the MHD runs, revealing that Hall-MHD dynamos 
can be more efficient (in the sense that they generate 
more magnetic energy and dissipate less total energy). 

Figure 6 shows the maximum value attained by the 
magnetic energy as a function of e for several simula- 
tions with different Reynolds numbers and scale sep- 
arations. When R = 560 the final amphtude reached 
by the magnetic energy is unknown. Given the strin- 
gent quadratic Courant-Friedrich-Lewy (CFL) condi- 
tion imposed by dispersive waves in Hall-MHD, sim- 
ulations were only carried up to saturation of the dy- 
namo. The maximum value of the energy in Figure 6 
is normahzed with the value obtained in an MHD run 
with the same Reynolds numbers and initial kinetic en- 
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Fig. 5. — Magnetic (below) and kinetic energy (above) 
as a function of time for two runs with e = 0.2 and 
(R = 300). 



ergy and helicity. As previously mentioned, consider- 
ing the three different regimes of the Hall dynamo dis- 
cussed in Mininni, Gomez, & Mahajan (2003b), we 
focus on the Hall-enhanced case. As the Reynolds 
numbers are increased, the efficiency of the Hall-MHD 
dynamo grows. Also, the value of e at which maxi- 
mum efficiency is obtained, decreases as the Reynolds 
numbers are increased. Note that the growth of the 
efficiency of the Hall-MHD dynamo with increasing 
scale separation was predicted analytically by Mininni, 
Gomez, & Mahajan (2002). The shift of most effi- 
cient e to smaller values as R increases is also obtained 
from analytical estimates (Mininni, Gomez, & Maha- 
jan 2004). 

5.2. Energy spectrum 

Figure 7 shows the kinetic, magnetic, and total 
energy spectra at different times for R = 300 and 
e = 0.1. Barring e, all the parameters and initial con- 
ditions in this run are the same as those corresponding 
to Figure 2. Therefore, a direct comparison between 
the evolution of both spectra can be made. During 
the first few time steps, the evolution is similar to the 
MHD run, with the entire magnetic spectrum growing 
at almost the same rate. The difference observed in 
Mininni, Gomez, & Mahajan (2003b), that the large- 
scale magnetic field is slightly larger than in its MHD 
counterpart, is now increased as a result of the larger 
Reynolds number and larger scale separation. 

The rate of increase of the large scale magnetic field 
changes in the presence of the Hall effect. While in 
MHD the build-up of this field proceeds on a resis- 
tive timescale, in Hall-MHD it grows faster. Note that 
the magnetic energy in the shell fc = 1 in Hall-MHD 
(e = 0.1) simulation att = 14.1 is a factor of 2 larger 
than the magnetic energy in the same shell in the MHD 
run att = 18.4. Also, the kinetic energy in the same 
shell is larger. This can be also observed in Figure 8, 
which shows the energy contained in the large scale 
magnetic field (in the shell k = 1) as a function of 
time for several values of e. However, the filling fac- 
tor 

((^B^^/(B2))1/2 is smaller as e is increased. At 

f = 14 the filhng factor is w 0.23 for e = 0, « 0.17 
for e = 0.1, and « 0.1 for e = 0.2, 

Note also that, while the MHD spectrum shows 
super-equipartition at small scales (the magnetic en- 
ergy is larger than the kinetic energy at large wave 
numbers), the Hall-MHD leads to equipartition at 
these scales. Both the evolution of the large scale 
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Fig. 6. — E,r,/E„i{e = 0) as a function of e for runs 
with R = 100 and R = 300. The values of cor- 
respond to the magnetic energy at the saturation level, 
while Em{e = 0) is the magnetic energy for an MHD 
run 




Fig. 7. — Mean kinetic energy spectrum (thick line), 
total energy spectrum (thick dashed line), and mag- 
netic energy spectrum at different times (e = 0.1 and 
R = 300). 



and small scale magnetic fields are therefore clearly 
affected by the Hall effect, even though the Hall effect 
operates effectively only at small scales. 

Figure 9 shows the spectrum for R = 300 and e = 
0.2. In this case aU wave numbers larger than kHaii = 
5 are affected by the Hall effect. However, the total 
energy spectrum in the saturated state seems to obey a 
Kolmogorov type law, although the Hall length scale 
is placed in the middle of the inertial range. 

Figure 10 shows the compensated energy spectrum 
E{k)/{e^/'^k^^^^) for higher spatial resolution runs 
with R = 560 and e = and 0.1, using 256"^ grid 
points. If the spectrum obeys a Kolmogorov type law, 
the compensated spectrum should be flat over a cer- 
tain range, and the amplitude of the spectrum in this 
range gives the Kolmogorov constant Ck- The mild 
hump observed for the total energy before entering the 
dissipative range, might be indicative of the presence 
of a "bottleneck effect" for the energy cascade, as was 
recently discussed by Haugen, Brandenburg & Dobler 
(2003). 

The spectra displayed in Figure 10 correspond to 
the era when the dynamo is saturated. The first exam- 
ple (Fig. 10(a)) represents the MHD limit (i.e. e = 0), 
while the second one (Fig. 10(b)) has kHaii = 10 
(e = 0.1). All the length scales smaller than l/knaii 
are expected to be dominated by the Hall effect. For 
the MHD spectrum, the Kolmogorov constant Ck ~ 
1.39, a value slightly larger than the one found by 
Haugen, Brandenburg & Dobler (2003), Ck « 1.3 
(see also Haugen, Brandenburg & Dobler (2004)), but 
smaller than the one obtained by Kida, Yanase, & 
Mizushima (1991) (Ck « 2.1). No clear change 
in the slope can be identified in the Hall-MHD case. 




Fig. 8. — Magnetic energy in the shell A; = 1 in Fourier 
space as a function of e and time. 
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and the spectrum is compatible with a Kolmogorov 
type law in its inertial range. However, the spectrum 
could really be a little shallower. Higher spatial reso- 
lution simulations are needed to settle this point. In 
this case the Kolmogorov's constant turns out to be 
Ck ~ 1-66, somewhat larger than its MHD counter- 
part. The Kolmogorov's dissipation wave numbers in 



70 [where 



the last stages of these runs are ki, 

k, = {{J') 

5.3. Magnetic helicity 

In Mininni, Gomez, & Mahajan (2003b), the Hall 
effect was observed to inhibit the creation of net mag- 
netic helicity by the dynamo process. This effect is en- 
hanced as we increase the Reynolds numbers. While 
the MHD dynamo is an efficient generator of magnetic 
helicity with most of this helicity concentrated in the 
larger scales, the Hall dynamo is somewhat sluggish; 
the growth of net magnetic helicity is slower and in 
some cases oscillates around zero (see Figure 11.) This 
result is in good agreement with theoretical estimates 
suggesting that in the presence of the Hall effect re- 
connection events are faster (Priest & Forbes 1998; Ji 
1999), and therefore dissipate less magnetic helicity 
(see Section 6). 

5.4. Kinetic helicity 

In MHD, relative kinetic helicity is known to 
change only slightly during dynamo action (Branden- 
burg 2001). In all these simulations, kinetic energy 
and kinetic helicity are injected at the same length 
scale by the stirring force acting at fc/orce = 3. For 
homogeneous hydrodynamic turbulence, kinetic helic- 
ity directly cascades to smaller scales. On dimensional 
grounds (Moffat 1978), the spectrum of kinetic helic- 
ity in the inertial range is also expected to follow Kol- 
mogorov's law (Cheng et al. 2003; Gomez & Mininni 
2004) 

Hk{k) = CHkheie^/^'k-^/^ , (21) 

where Ch is another Kolmogorov constant, and khei 
is the scale where kinetic helicity is injected. Using 
the inertial range kinetic energy expression {Ck is the 
standard Kolmogorov's constant ) 



Ek{k) 



(22) 



the ratio between the total kinetic helicity and kinetic 
energy in hydrodynamic turbulence come out to be 



Ek 



JU-u>dV 

JU^dV 



A 






k"='' t=0 




t=2.9 : 




. . t=3.9 : 




^\ 1 = 5,4 








































— - 















Fig. 9. — Mean kinetic energy spectrum (thick line), 
total energy spectrum (thick dashed line), and mag- 
netic energy spectrum at different times (e = 0.2 and 
R = 300). 





Fig. 10. — Total compensated energy spectrum 
^5/3^-2/3^^ (thick line), kinetic energy spectrum 
(thin line), and magnetic energy spectrum (dashed 
line) at t = 5.5 for a run with (a) e = 0, and (b) e = 0. 1 
{R = 560, 256^ grid points). 



khel 



(23) 
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where khei ~ k force in our case. What happens after 
we introduce the magnetic seed? 

Figure 12 shows the ratio Hi~ / for the MHD and 
the Hall-MHD runs with R = 300. The evolution of 
the relative kinetic helicity is similar. In the MHD run 
Hk/Ek « 2.5, a value close to k force = 3. Note that, 
although the total kinetic energy (and the kinetic heUc- 
ity) decreases during the time evolution as a result of 
the increasing Lorentz force, the ratio Hk / Ek remains 
nearly constant. On the other hand, this ratio grows 
with e in the HaU-MHD runs. The growth takes place 
just after the exponentially growing stage, when a large 
scale magnetic field is developing in the box. This re- 
sult suggests that a new source of kinetic helicity has 
appeared, and is in good agreement with theoretical es- 
timates suggesting that the Hall effect introduces hand- 
edness in the fluid motions (Mininni, Gomez, & Ma- 
hajan 2004). However, we want to point out that this 
handedness does not, by itself, generate a net a-effect. 

As previously mentioned, the relative kinetic he- 
Ucity Hk/{{U^) (a;2)) in Hall-MHD also shows the 
same behavior, changing from a relative kinetic helic- 
ity of 0.4 when e = up to about 0.6 when e = 0.1. 

6. LARGE SCALE MAGNETIC FIELD GEN- 
ERATION 

In this section we present MHD and Hall-MHD re- 
sults for an external force located at k force = 10. Sim- 
ulations were carried with v = r] = 0.02 and 128^ grid 
points. In this case the Reynolds numbers are smaller 
(R f» 220 and R\ « 20), and the turbulence is weaker, 
since there are not enough modes in Fourier space for a 
direct cascade to develop properly. On the other hand, 
there is more room for the generation of a large scale 



magnetic field through inverse cascade. Simulations 
were performed for e = 0, 0.1, and 0.2, correspond- 
ing respectively to the MHD case, kHaii = 10 and 
knaii = 5. 

Figure 13 shows the evolution of the magnetic and 
kinetic energy in the MHD simulation. After satura- 
tion, a large scale magnetic field grows on a resistive 
time scale, as will be shown in the energy spectrum and 
was previously observed in large scale dynamo simu- 
lations (Brandenburg 2001). Note that the system fi- 
nally reaches a state of super-equipartition, i.e. a level 
of magnetic energy which is larger than the kinetic en- 
ergy. Figure 14 shows its counterpart when e = 0.1 
and 0.2. In the simulation with e = 0.1, the growth of 
magnetic energy after the saturation is clearly faster, 
and the system reaches a final state with more mag- 
netic energy than in the MHD run. 

Figures 15 and 16 show the energy spectrum at dif- 
ferent times for runs with e = and e = 0.1, respec- 
tively. Note that in both simulations, the magnetic en- 
ergy at intermediate scales (2 < fc < 8) starts to de- 
cay after saturation (t > 3.5), while magnetic energy 
at the largest scale (k = 1) keeps growing. As a re- 
sult, after < = 8 the system reaches the state of super- 
equipartition in the MHD case. This is even more clear 
in Fourier space (Figure 15), where magnetic energy in 
the shell A: = 1 is two orders of magnitude larger than 
kinetic energy at t = 17.5. An excess of magnetic 
energy can also be observed at small scales. 

In the Hall-MHD case, we observe again a faster 
growth of the large scale magnetic field (Figure 16), 
but with a final state with super-equipartition only 
at large scales. Moreover, the kinetic energy in the 
sheU fc = 1 is one order of magnitude larger than 
in the MHD case, and at smaU scales we obtain sub- 



1.0 




0.5 




0.0 






i" -0.5 


r F=n \ \ 


-1.0 


E = 0.066 

: e = 0.1 : 


-1.5 


- .... 6 = 0.2 - 


-2.0 





2 4 6 8 10 12 

Time 

Fig. 11. — Magnetic helicity for e = 0, 0.066, 0.1, and 
0.2 (R = 300). 




2.2 - 

2.0 E , , , , I , , , , ^ 

5 10 

Time 



Fig. 12.— Hk/Ek for e = 0, 0.066, and 0.1. 
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Fig. 13. — Magnetic (below) and kinetic energy 
(above) as a function of time for e = and k force = 
10. 



equipartition between magnetic and kinetic energy. 

Figure 17 shows the magnetic helicity spectrum 
as a function of time for different values of e. As 
mentioned in Section 4, the a effect creates magnetic 
helicity of a sign opposite to that of kinetic helicity 
at large scales. This effect is balanced by the cre- 
ation of an opposite amount of magnetic helicity at 
small scales. Therefore, the diffusion preferentially 
destroys the short-scale magnetic helicity in reconnec- 
tion events, leaving a net helicity of opposite sign at 
large scales (Brandenburg 2001). Mininni, Gomez, 
& Mahajan (2003b) suggested that while the Hall- 
MHD dynamo process also creates equal and oppo- 
site amounts of magnetic hehcity at large and at small 
scales, the dissipation of magnetic helicity at small 
scales is less efficient as e is increased. Figure 17 
shows that when the Hall effect is present, even at late 
times an excess of positive magnetic helicity at small 
scales (A; > 10) can be readily identified in the spectra. 

7. DISCUSSION 
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Fig. 14. — Magnetic (below) and kinetic energy 
(above) as a function of time for two runs with e = 0.2, 
O.l and Oik force = 10). 



In this paper we have presented the results of di- 
rect numerical simulations of turbulent dynamo action 
in Hall-MHD. We find that with increasing Reynolds 
number and scale separation, the Hall MHD dynamo 
works more efficiently when the Hall length is close 
but larger than the dissipation scale (Hall-enhanced 
regime). For larger values of e the Hall MHD dynamo 
is less efficient. In addition, the value of e (which mea- 
sures the strength of the Hall term) at which the dy- 
namo is most efficient decreases at higher Reynolds 
numbers. 

An acceleration of the process responsible for the 
growth of a large scale magnetic field is observed at 
moderate values of e. Although these simulations are 
made at Reynolds numbers which are far away from 
reahstic values for astrophysical plasmas, the results 
obtained are encouraging; the dynamos tend to work 
better at high Reynolds numbers 

By calculating the magnitude and nature of the gen- 
erated magnetic field as the ampUtude of the Hall term 
is varied, we obtain new evidence showing that the 
Hall dynamo can be fundamentally different from its 
classical MHD counterpart. 
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total energy spectrum (thick dashed hne), and mag- 
netic energy spectrum at different times (e = and 

kforce — 10)- 




0.001 , , , , , . . , ■ ' 

1 10 

k 

Fig. 16. — Mean kinetic energy spectrum (thick Une), 
total energy spectrum (thick dashed Une), and mag- 
netic energy spectrum at different times (e = 0.1 and 

kforce — 10)- 



reviewing the manuscript, and to the Abdus Salam In- 
ternational Centre for Theoretical Physics, were the 
initial stages of this study were performed. Research 
of SMM was supported by US DOE contract DE- 
FG03-96ER-54366. Research of DOG and PDM has 
been partially funded by grant X209/01 from the Uni- 
versity of Buenos Aires and by grant PICT 03-9483 
from ANPCyT. PDM is a fellow of CONICET, and 
DOG is a member of the Carrera del Investigador Ci- 
entifico of CONICET. 

REFERENCES 

Balbus, S.A. & Terquem, C. 2001, ApJ, 552, 235 

Blackman, E.G. & Field, G.B. 1999, ApJ, 521, 597 

Braginskii, S.l. 1965, Rev. Plasma Phys., 1, 205 

Brandenburg, A., & Subramanian, K. 2000, Astron. & 
Astrophys., 361, L33 

Brandenburg, A. 2001, ApJ, 550, 824 

Canuto, C, Hussaini, M.Y., Quarteroni, A., & Zang, 
TA. 1988, Spectral Methods in Fluid Dynamics 
(BerUn: Springer- Verlag) 

Galanti, B., Kleeorin, N., & Rogachevskii, I. 1994, 
Phys. of Plasmas, 1,3843 

Cattaneo, F. & Hughes, D.W. 1996, Phys. Rev. E, 54, 
R4532 

Chen, Q., Chen, S., & Eyink, G. 2003, Phys. Fluids, 
15,361 

Gomez, D. & Miniimi, P. 2004, Phys. A, in press. 

Haugen, N.E.L., Brandenburg, A., & Dobler, W. 2003, 
ApJ, 597, L141 

Haugen, N.E.L., Brandenburg, A., & Dobler, W. 2004, 
Phys. Rev. E, 70,016308 

Heintzmaim, H. 1983, JETP, 57, 251 

Helmis, G. 1968, Mber. Dtsch. Akad. Wiss. Berlin, 10, 

280 

Helnois, G. 1971, Beitr. Plasma Physik, 1 1 417 

Ji, H. 1999, Phys. Rev. Lett, 83, 3198 

Kida, S., Yanase, S., & Mizushima J. 1991, Phys. Flu- 
ids, A3, 457 



12 



Krause, F. & Radler, K.-H., 1980, Mean-Field Mag- 
netohydrodynamics and Dynamo Theory, (GDR: 
Pergamon Press) 

Meneguzzi, M., Frisch, U., & Pouquet, A. 1981, 
Phys. Rev. Lett., 47, 1060 

Miniimi, PD., Gomez, D.O., & Mahajan, S.M. 2002, 
ApJ, 567, L81 

Mininni, RD., Gomez, D.O., & Mahajan, S.M. 2003a, 
ApJ, 584, 1120 

Mininni, RD., Gomez, D.O., & Mahajan, S.M. 2003b, 
ApJ, 587, 472 

Mininni, RD., Gomez, D.O., & Mahajan, S.M. 2004 
(in preparation). 

Mimov, V.V., Hegna, C.C., & Prager, S.C. 2003, 
Plasma Phys. Rep., 29, 566 

Moifat, H.K. 1978, Magnetic field generation in elec- 
trically conducting fluids (Cambridge: Cambridge 
University Press) 

Muslimov, A.G. 1994, MNRAS, 267, 523 

Pope, S.B. 2000, Turbulent flows (Cambridge: Cam- 
bridge University Press) 

Pouquet, A., Frisch, U., & Leorat, J. 1976, J. Fluid 
Mech., 77, 321 

Priest, E. & Forbes, T. 1998, Magnetic Reconnection 
(Cambridge: Cambridge Univ. Press) 

Rheinhardt, M. & Geppert, U. 2002, Phys. Rev. Lett., 

88, 101103 

Sano, T. & Stone, J.M. 2002, ApJ, 570, 314. 
Zweibel, E.G. 1988, ApJ, 329, 384 





Fig. 17. — Magnetic helicity spectra for: (a) e = 0, 
and (b) e = 0.2 (R = 300). 



This 2-column preprint was prepared witli tlie AAS lATj5C macros 
v5.0. 



13 



